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ABSTRACT 

Optimization of rotorcraft flowfields using an adjoint method generally requires a time-dependent implementation of the 
equations. The current study examines an intermediate approach in which a subset of rotor flowfields are cast as steady 
problems in a noninertial reference frame. This technique permits the use of an existing steady-state adjoint formulation with 
minor modifications to perform sensitivity analyses. The formulation is valid for isolated rigid rotors in hover or where the 
freestream velocity is aligned with the axis of rotation. Discrete consistency of the implementation is demonstrated using 
comparisons with a complex-variable technique, and a number of single- and multi-point optimizations for the rotorcraft 
figure of merit function are shown for varying blade collective angles. Design trends are shown to remain consistent as the 
grid is refined. 


NOTATION 

C Aerodynamic coefficient 

Cq Rotor torque coefficient 

C T Rotor thrust coefficient 

D Vector of design variables 

E Total energy per unit volume, modulus of elasticity 
F / ,F V In viscid and viscous flux tensors 

FM Rotorcraft figure of merit function 
I Identity tensor 

K Elasticity coefficient matrix 

L Lagrangian function 

N Number of composite objective functions 

Q Vector of conserved variables 

R Spatial residual vector 

S Source term vector 

S Control volume surface area 

T Temperature 

V Volume of control volume 
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X Vector of grid coordinates 

/ Objective function 

g Real-valued function 

h Step size 

i A 

i, j , k Indices 
k Thermal conductivity 

m Number of constraint function components 

n Outward-pointing normal vector 

n Number of objective function components 

p Pressure, exponent 

r Position vector 

t Time 

iij Cartesian directional displacements 

u,v,w Cartesian components of velocity 
x Independent variable 

Xj Cartesian coordinate directions 

A j' Flowfield adjoint variable 

A g Grid adjoint variable 

0 Blade collective setting 

12 Angular velocity vector 

a , co Weights 
ti Strain tensor 

1 ) Spanwise station 



A , p Lame constants 
p Density, K-S multiplier 

r Viscous stress tensor 

V Poisson’s ratio 
mp Denotes multi-point quantity 

surf Denotes surface quantity 

* Denotes target quantity, optimal quantity 

INTRODUCTION 

Application of high-fidelity computational fluid 
dynamics (CFD) has become commonplace in the fixed- 
wing aerospace community. Software packages that solve 
the Euler and Reynolds-averaged Navier-Stokes equations 
on both structured and unstructured grids are now used 
routinely in the analysis and design of new configurations. 
Moreover, as algorithms and computer hardware have 
continued to mature, the use of formal design optimization 
techniques coupled with CFD methods has become viable 
for large-scale problems in aerospace design. 

The application of high-fidelity CFD tools to the 
analysis and design of full rotorcraft configurations is 
considerably more challenging. Such flowfields are 
inherently unsteady, frequently involve fluid velocities 
ranging from quiescent to transonic flow, and typically 
require the simulation of complex aerodynamic and 
aerostructural interactions between dynamic vehicle 
components. Recent literature suggests the use of high- 
fidelity CFD methods in this regime is growing, but the 
computational cost required to capture the necessary spatial 
and temporal scales of a typical rotorcraft flowfield remains 
considerable. [1H91 

In the field of gradient-based design, adjoint methods 
are known to provide an extremely efficient means for 
computing sensitivity information. The cost of such methods 
is equivalent to the expense associated with solving the 
analysis problem and is independent of the number of design 
variables. Adjoint methods can also be used to perform 
mathematically rigorous mesh adaptation and error 
estimation. Significant success has been reported for the 
application of these techniques to steady problems; for 
example, see Refs. [ 10]-[ 1 3] and the efforts cited in Ref. 
[14]. 

In general, optimization and mesh adaptation for large- 
scale rotorcraft flows using adjoint methods require a time- 
dependent implementation of the equations. Considerable 
effort by a number of research groups is being focused in 
this area, and examples of the use of such approaches have 
just recently emerged. [151[171 Despite the algorithmic 
efficiency however, the computational cost of these general 
time-dependent approaches can be considerable, and the 
application of such methods to practical problems of 
engineering interest may remain prohibitively expensive for 
some time. 

The goal of the current work is to develop, implement, 
and demonstrate an adjoint-based design capability for rotor 
configurations for which the analysis problem may be cast as 
a steady problem in a noninertial reference frame. This 
approach permits the use of an existing steady-state adjoint 


formulation with minor modifications to perform sensitivity 
analyses. The resulting formulation is valid for isolated rigid 
rotors in hover or where the freestream velocity is aligned 
with the axis of rotation. 


FLOW EQUATIONS 


The governing equations for the flowfield are the 
compressible, perfect gas Reynolds-averaged Navier-Stokes 
equations written in a reference frame rotating with a 
constant angular velocity O : 

^ + d) (F ; - -¥ v )-ndS = S , (1) 

at J gy 

where Q is the vector of volume-averaged conserved 
variables Q = [p,pu,pv,pw, £] , n is an outward pointing 
unit normal, and V is the control volume bounded by the 
surface dV . The inviscid and viscous flux tensors are given 
by 
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The source term S represents a Coriolis effect due to the 
rotating frame of reference: 

r o i 
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Here, u is the absolute velocity vector u = [u,v, vv] , r is the 
position vector relative to the axis of rotation, and r is the 
viscous stress tensor. The equations are closed with the 
perfect gas equation of state and an appropriate turbulence 
model for the eddy viscosity. For rotorcraft simulations, the 
formulation described here is applicable to rigid rotor 
geometries in either a hover condition or 

ascending/descending flight, where the freestream velocity 
vector is parallel to the angular velocity vector . 

References [18]-[21] describe the flow solver used in 
the current work. The code can be used to perform 
aerodynamic simulations across the speed range and an 
extensive list of options and solution mechanisms is 
available for spatial and temporal discretizations on general 
static or dynamic mixed-element unstructured meshes which 
may or may not contain overset grid topologies. 

In the current study, the spatial discretization uses a 
finite-volume approach in which the dependent variables are 
stored at the vertices of single-block tetrahedral meshes. 
Inviscid fluxes at cell interfaces are computed using the 
upwind scheme of Roe, [221 and viscous fluxes are formed 
using an approach equivalent to a central-difference 
Galerkin procedure. The eddy viscosity is modeled using the 
one-equation approach of Spalart and Allmaras [231 with the 
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source term modification proposed by Dacles-Mariani. [241 
For the steady-state flows (relative to the noninertial 
reference frame) described in this study, temporal 
discretization is performed using a backward-Euler scheme 
with local time stepping. Scalable parallelization is achieved 
through domain decomposition and message passing 
communication. 

An approximate solution of the linear system of 
equations formed within each time step is obtained through 
several iterations of a multicolor Gauss-Seidel point-iterative 
scheme. The turbulence model is integrated all the way to 
the wall without the use of wall functions. The turbulence 
model is solved separately from the mean flow equations at 
each time step with a time integration and linear system 
solution scheme identical to that employed for the mean 
flow equations. 


GRID EQUATIONS 


To deform the interior of the computational mesh as the 
surface grid evolves during a shape optimization procedure, 
the mesh is assumed to obey the linear elasticity equations of 
solid mechanics. These relations can be written as 


/ 3 

d) A V T-ntfS + d) 2/jx ■ ndS = 0 , 

J dv J dv 
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is the strain tensor, m ; - is the displacement vector in each of 
the Cartesian coordinate directions x i , and A and fi are 
material properties of the elastic medium. The quantities A 
and fi are related to Young’s modulus E and Poisson’s 
ratio V through the following: 


A = 


vE 

(\ + V)(\- 2 v) 


(7) 


and 


E 

2(1 + v ) ' 


( 8 ) 


The system is closed with the specification of two of the four 
parameters A , // , E , and V . In the current 

implementation, E is taken as inversely proportional to the 
distance from the nearest solid boundary, while Poisson's 
ratio is taken uniformly as zero. This approach forces cells 
near boundaries to move in a nearly rigid fashion, while cells 
far from the boundaries are allowed to deform more freely. 
The system of equations is solved using GMRES [251 with 
either a point-implicit or ILU(O) preconditioning technique 
as described in Refs. [21] and [26]. 


DISCRETE ADJOINT EQUATIONS 

To derive the discrete adjoint equations, it is useful to 
introduce a compact notation for the governing equations 
outlined above. The spatial residual vector R of Eq. (1) is 
defined as 


R = d) (F ; - - F v ) • hdS - S . (9) 

J dv 

Furthermore, the linear system of equations given by Eq. (5) 
can be written as 

KX = X sulf , (10) 

where K is the elasticity coefficient matrix resulting from 
the discretization of Eq. (5), X is the vector of grid point 
coordinates, and X su/ y is the vector of known surface grid 
point coordinates, complemented by zeros for all interior 
coordinates. 

Following the approach taken in Ref. [11], a Lagrangian 
function can be defined as follows: 

L(D,Q,X,A / ,A g ) = /(D,Q,X) 


+A jR(D,Q, X) + A^ (KX — X iH) y ) , 


( 11 ) 

where D represents a vector of design variables, / is an 
objective function, and Ay and A ? are adjoint variables 
multiplying the residuals of the flow and grid equations. In 
this manner, the governing equations may be viewed as 
constraints. 

Differentiating Eq. (11) with respect to D and equating 
the 5Q/5D and 3X/3D coefficients to zero leads to the 
discrete adjoint equations for the flowfield and grid, 
respectively: 
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(13) 


The remainder of the terms in the linearized Lagrangian can 
be grouped to form an expression for the final sensitivity 
vector: 


dL _ df j 3R j dX sur f 

dD 3D J 3D g 3D 


(14) 


Eqs. (12) and (13) provide a very efficient means for 
determining discretely consistent sensitivity information. 
The expense associated with solving these equations is 
independent of D , and is similar to that of the solution of 
the governing equations. Once the solutions for A y and 
A g have been determined, the desired sensitivities may be 
calculated using Eq. (14), whose computational cost is 
negligible. 

A discrete adjoint implementation has been developed 
in Refs. [11], [17], [20], [26], and [27] for the flow solution 
method described above. The flowfield adjoint equations are 
solved in an exact dual fashion which ultimately guarantees 
an asymptotic convergence rate identical to the primal 
problem and costate variables which are discretely adjoint at 
every iteration of the solution process. The grid adjoint 
equations are solved using GMRES in a manner identical to 
that done for Eq. (5). To accommodate the noninertial 
reference frame used in the current study, minor 
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modifications have been made to include the effects of the 
mesh speeds and Coriolis terms. 

DESIGN METHODOLOGY 
Design Variables 

The implementation described in Ref. [11] is 
sufficiently general that the user is able to employ a 
geometric parameterization scheme of choice, provided the 
associated linearizations required by the adjoint method 
described above are also available. For the current study, the 
grid parameterization scheme described in Ref. [28] is used. 
This approach can be used to define very general shape 
parameterizations of existing grids using a set of aircraft- 
centric design variables such as camber, thickness, shear, 
twist, and planform parameters at various locations on the 
geometry. The user also has the freedom to directly associate 
two or more design variables to create more general 
parameters. In the current work, this option is used to link 
several twist variables across the span of a rotor blade to 
create a single twist variable that is used to prescribe the 
blade collective setting 0 . In the event that multiple bodies 
of the same shape are to be designed - as in the case of rotor 
geometries - the implementation allows a single set of 
design variables to be used to simultaneously define such 
bodies. In this fashion, the geometry of each body remains 
consistent throughout the course of the design. 

Objective and Constraint Functions 

The implementation described in Ref. [11] permits 
multiple objective functions f- t and explicit constraints Cj of 
the following form, each containing a summation of and 
nij individual components, respectively: 

fi = ioJ k (C k -C k ^ (15) 

k = 1 

and 

c j = Ta> k (C k -C* k ) p * . (16) 

k=\ 

Here, 0 \ represents a user-defined weighting factor, C k is 
an aerodynamic coefficient such as total drag or the pressure 
or viscous contributions to such quantities, and p k is a user- 
defined exponent. The (*) superscript indicates a user- 
defined target value of C k . Furthermore, the user may 
specify which boundaries in the grid to which each 
component function applies. 

Design Points and Optimization Strategies 

The current implementation supports an arbitrary 
number of user-specified design points where objective and 
constraint functions may be posed. Each design point may 
be characterized by a variation of basic flowfield quantities 
such as the Mach number, or a more general characteristic 
such as the computational grid appropriate for each 


individual design point. In the current study, each blade 
collective setting 0 requires a different grid and therefore 
represents a different design point. 

To perform multi-point optimization, three methods are 
considered. The first two approaches are unconstrained 
formulations where individual objective functions /) are 
posed at each design point, from which an overall composite 
objective function f mp is constructed. The third approach is 
a constrained formulation. 

The first method used to form the composite objective 
function f mp defines a linear combination of /'- : 

fmp =a \f\ +a 2fl +a lh + --- + a NfN' (l 7 ) 


where N is the total number of design points and CCj is a 
constant weighting factor applied to each individual f; . In 
the current study, all values of a t are chosen to be 1.0. 

The second approach used to define f mp is based on the 
technique described in Refs. [29] and [30]. In this approach, 
the objective functions /'- from each design point are 
combined using the Kreisselmeier-Steinhauser function to 


form f, 


mp • 


■fmp /max ^ 


P 


N 

£■ 

. 1=1 


cax) 


( 18 ) 


The quantity / max is defined as the maximum value over all 
fi and the value p is a user-defined constant taken to be 
20.0. 1301 Although not considered here, this approach also 
has the added benefit of being able to convert constrained 
optimization problems into unconstrained problems by 
including explicit constraints in the formulation of Eq. (18). 

The third multi-point formulation considered is based 
on a constrained formulation. In this approach, the objective 
function to be minimized is defined at a single design point, 
while the objective functions defined at the other design 
points are instead treated as explicit constraints on the 
optimization problem. 

The multi-point approaches used here are common in 
obtaining point solutions to multi-objective optimization 
problems via scalarization of the multiple objectives. The 
difficulty is that out of the range of many possible solutions 
only one is obtained by setting some parameters heuristically 
and externally, e.g., the weights of the composite scalar 
objective. Since the current focus is the interaction of adjoint 
methods with design optimization, in principle, these simple 
strategies adopted here suffice, but it is noted that the related 
areas of robust and multi-objective design are extensive and 
active. The investigation of more sophisticated optimization 
strategies is relegated to future work. 

For unconstrained problems, the optimization package 
described in Ref. [31] is used to minimize the specified 
objective function. In these cases, the optimizer is allowed to 
perform up to 20 design cycles or 30 function evaluations, 
whichever occurs first. The optimization algorithm considers 
the design converged and exits if it believes the following 
stopping tolerance is met: 
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f f <lxlO~ 5 , (19) 

|/| 

where / is the objective function value at the optimal 
solution to the design problem and is not known a priori. 

The package outlined in Ref. [32] is used for problems 
where explicit constraints are present. The optimization 
algorithm is allowed to perform a maximum of 20 design 
cycles, and considers the design converged and exits if it 
believes the current objective function matches the value at 
the optimal solution to four significant digits. Constraints are 
considered satisfied if their values do not exceed the 
specified bounds by 0.5% of the bound value. The design at 
the initial choice of D is not required to satisfy the 
constraints; if needed, the optimizer will attempt to locate a 
feasible starting point on its own. 

TEST CASE 

Demonstration optimizations are computed using the 
three-bladed Tilt Rotor Aeroacoustics Model (TRAM) 
described in Refs. [33] and [34] and shown in Fig. 1. The 
optimizations are performed for a hover condition 
corresponding to collective settings 0 = 10°, 0 = 12°, and 
0 = 14° . The tip Mach number is 0.62 and the Reynolds 
number is 2.1 million based on the blade tip chord. The 
mesh used for the design studies contains 5,048,727 nodes 
and 29,802,252 tetrahedral elements and is designed for the 
0 = 14° setting. Grids for the 0 = 10° and 0 = 12° settings 
are obtained through elastic deformations of the baseline 
mesh. The surface grid for one of the blades is shown in 
Figure 2. All of the grids have been generated using the 
approach outlined in Ref. [35]. 

A geometric parameterization has been developed for 
the baseline blade geometry as shown in Fig. 3. The 
approach yields a total of 44 active design variables 
including 20 variables to control the blade thickness and 24 
variables to control the blade camber. The root section of 
each blade is held fixed. Bounds on the design variables 
have been initially chosen with the intent of preventing non- 
physical surface shapes; further constraints on the minimum 
thickness will be described in a later section. The 
parameterization also allows for blade planform variations as 
well as local twist and shearing deformations; however, 
these are held fixed in the current study. 

For each design point, a single objective or constraint 
function is used, where co = 1 , p = 2 , and C is defined as 
the square of the commonly-used rotorcraft figure of merit 
function, composed of the rotor thrust and torque 
coefficients: 

C 3 

C = FM 2 = — . (20) 

K 

The square has been introduced to avoid the possibility of a 
square root of a negative thrust value appearing in the 
linearized form of the objective function. In all cases, the 
value of C is chosen to be 2.0, which is considerably larger 


than the baseline value at each of the collective settings 
examined here, as well as the theoretical maximum value of 

1 . 0 . 

All computations have been performed using 75 3.0 
GHz dual-core Pentium IV processors with gigabit ethernet 
connections. A typical design cycle requires a single 
function and gradient evaluation for the current value of D . 
A function evaluation in this context consists of an 
evaluation of the surface parameterization for each blade, a 
solution of Eq. (5) to deform the interior of the mesh 
according to the current surface grid, and a solution of the 
flow equations, Eq. (1). Using the adjoint approach outlined 
above, a gradient evaluation requires a solution of the 
flowfield adjoint equations, Eq. (12); a solution of the mesh 
adjoint equations, Eq. (13); an evaluation of the linearized 
surface parameterization for each blade; and finally, an 
evaluation of the gradient expression given by Eq. (14). This 
combined procedure for obtaining a single function and 
gradient vector for a given collective setting 0 takes 
approximately 2.5 wallclock hours using the stated 
hardware. The convergence criteria used for each of the 
solvers has a direct impact on this efficiency. Finally, the 
time required to solve Eqs. (1) and (12) tends to decrease 
towards the end of an optimization as the design converges 
and solution restarts become more effective. 


ACCURACY OF IMPLEMENTATION 


To verify that a discretely consistent implementation of 
Eqs. (12)-( 14) has been achieved, results are compared with 
those obtained using an independent approach based on the 
use of complex variables. This technique was originally 
suggested in Refs. [36] and [37], and was first applied to a 
Navier-Stokes solver in Ref. [38]. In this approach, a Taylor 
series with a complex step size ih is used to derive an 
expression for the first derivative of a real-valued function 

g(x): 


g\x) 


Im[g(jc + i/i)] 
h 



( 21 ) 


The primary advantage of this approach is that true second- 
order accuracy may be obtained by selecting step sizes 
without concern for subtractive cancellation error typically 
present in real-valued divided differences. This capability 
can be immediately recovered at any time for the baseline 
solvers used in this study through the use of an automated 
scripting procedure as outlined in Ref. [39]. 

A coarse mesh consisting of 144,924 nodes and 848,068 
tetrahedral elements is used to demonstrate the accuracy of 
the implementation for fully turbulent flow at the stated test 
conditions and a 0 = 14° collective setting. Sensitivity 
derivatives of the figure of merit with respect to several 
shape parameters located at the midspan location of each 
blade are computed using the discrete adjoint 
implementation. Results are compared with values obtained 
using the complex variable method, where a step size 

-5 A 

/z = lxl0 has been chosen. All equation sets are 
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converged to machine precision using 16 processors and 
results are shown in Table 1. The sensitivity derivatives 
computed using the two methods are in excellent agreement. 

RESULTS 

Single-Point Designs 

The first set of results is a single-point design at each of 
the chosen blade collective settings. The history for the 
figure of merit during the course of each design is shown in 
Fig. 4. For each collective setting, the figure of merit 
increases quickly during the early portion of the 
optimization, after which further gains are minimal. The 
initial and final figures of merit for each 0 are listed in 
Table 2. Improvements range from 4% to just under 8%, 
with smaller improvements at the higher collective settings. 
An expanded view of the resulting blade shape for each 
design is shown in Fig. 5, where the blades have each been 
rotated to the 0 = 14° setting for comparison purposes. The 
design changes are similar at each collective setting: the 
camber has been increased across the majority of the span, 
while the thickness has been reduced. Of particular interest 
is the blade trailing edge, where each design has reduced the 
blade thickness to a numerically valid but physically 
infeasible dimension. Where the thickness is fixed at the 
blade tip, the optimization has increased the camber for the 
0 = 10° setting, while decreasing it for the other two 
collectives, most notably for the 0 = 14° setting. 

Single-Point Designs with Thickness Constraints 

In an effort to achieve a more practical blade design in 
the trailing edge region, the previous set of test cases is 
repeated. However, constraints are now placed on the 
thickness variables to enforce the original blade thickness as 
a lower bound. The results using this approach are shown in 
Figs. 6 and 7. As before, the figure of merit for each 
collective setting is increased rapidly during the initial 
portion of the optimization. Examination of the blade cross- 
sections shows that the thickness of the baseline airfoil shape 
has been maintained as a lower bound. Differences between 
the designs at the various collective settings can be readily 
seen at the ?) = 0.40 station and the blade tip. Table 3 shows 
the figure of merit results for each collective setting. The 
improvements are less than those observed where blade 
thinning was allowed, ranging from almost 3% to 5.6%, with 
the largest improvements again taking place at the lower 
collective settings. 

Multi-Point Designs 

To evaluate the implementation for multi-point 
optimization problems, designs are performed using the 
three strategies outlined earlier. For the approach involving 
explicit constraints, the objective function is defined at the 
0 = 14° setting, while the functions defined at the other two 
collective settings serve as constraints. The lower bounds 


placed on these constraints correspond to minimum figures 
of merit of 0.71 and 0.73 at the 0 = 10° and 0 = 12° 
settings, respectively. These choices represent moderate 
increases over the baseline figure of merit at each 0 based 
on the single-point design results. Note that since these 
constraints are not satisfied by the initial blade geometries in 
this approach, the optimization procedure must locate the 
feasible region during the course of the design. The 
minimum thickness constraint is also enforced for each of 
the three multi-point approaches. 

The convergence history for the approaches based on 
Eqs. (17) and (18) are shown in Figs. 8 and 9, respectively. 
The two approaches yield comparable behavior for the 
figure of merit at each 0 . The final values given in Tables 4 
and 5 are also similar, although slightly higher for the 
approach based on the linear combination of individual 
objectives. 

The convergence for the constrained approach is shown 
in Fig. 10. The blade design satisfies the constraints at 
0 = 10° and 0 = 12° after the first design cycle, and the 
overall convergence for each collective setting is similar to 
the previous cases. However, it should be noted that for this 
particular case the optimization procedure was terminated 
early due to queue limitations on the computational 
platform. The procedure could be restarted if desired, but 
this has not been pursued here. Table 6 shows that the final 
blade design using this approach has figure of merit values 
that are comparable to the other multi-point approaches. 

Although the final figures of merit obtained through 
each of the multi-point methods are similar, the differences 
in the optimized blade geometries are striking, as shown in 
Fig. 11. An investigation of the off-design performance for 
each blade geometry and introduction of multidisciplinary 
interactions in the design process are logical next steps but 
beyond the scope of the current work. 

Effect of Grid Refinement 

A grid refinement study is performed using the initial 
and final geometries resulting from the multi-point 
optimizations described above. For these computations, a 
refined grid consisting of 12,662,080 nodes and 87,491,279 
tetrahedra has been constructed and parameterized in a 
manner consistent with the baseline grid. The final design 
variables established in the multi-point optimizations are 
applied to the refined grid and a single analysis is performed 
for each geometry to evaluate the resulting figure of merit. 
Results for the refined grid are included in parenthesis in 
Tables 4-6 beneath the values for the baseline grid. Although 
the magnitude of the design improvements varies slightly 
with grid density, the results on the refined grid show similar 
trends in all cases as compared with the baseline mesh. 

CONCLUDING REMARKS 

A discrete adjoint-based methodology for performing 
design optimization of isolated rotor problems which appear 


6 



as steady flows in a noninertial reference frame has been 
developed and implemented. The accuracy of the 
linearization has been established using comparisons with an 
independent approach based on the use of complex 
variables. A series of single- and multi-point designs at 
several blade collective settings showed improvements in the 
figure of merit function for both unconstrained and 
constrained problem formulations. Design trends were 
shown to remain consistent with grid refinement. 

Ongoing efforts are focused on a general time- 
dependent adjoint-based optimization capability for 
rotorcraft as well as other aerospace configurations 
characterized by unsteady flowfields. The efficiency of such 
an implementation should be compared with that of the 
present approach as well as other techniques such as time- 
periodic formulations. 
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Figure 1. Surface geometry for TRAM rotor. 



Figure 2. Typical blade surface grid used for design 
computations. 
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Figure 3. Design variable and radial blade locations. 



Figure 4. Figure of merit histories for optimizations with 
no thickness constraints included. 





Figure 5. Blade cross-sections at various radial stations 
before and after optimization with no thickness 
constraints included. The vertical scale has been 
exaggerated and all blades have been rotated to the 
0 = 14° collective setting for comparison. 



Figure 6. Figure of merit histories for optimizations with 
thickness constraints included. 



Figure 7. Blade cross-sections at various radial stations 
before and after optimization with thickness constraints 
included. The vertical scale has been exaggerated and all 
blades have been rotated to the 0 = 14° collective setting 
for comparison. 



Figure 8. Figure of merit histories for multi-point 
optimization based on the linear combination of objective 
functions given by Eq. (17). 
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Figure 9. Figure of merit histories for multi-point 
optimization based on the KS function given by Eq. (18). 
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Figure 10. Figure of merit histories for multi-point 
optimization based on the explicitly constrained 
approach. Arrows indicate feasible side of constraints. 



Figure 11. Blade cross-sections at various radial stations 
before and after multi-point optimization using the three 
different strategies considered. The vertical scale has 
been exaggerated and all blades have been rotated to the 
0 = 14 ° collective setting for comparison. 


Table 1. Comparison of figure of merit sensitivity 
derivatives obtained using adjoint and complex variable 
approaches. “A” denotes adjoint result, “C” denotes 
complex-variable result. 

Design Variable d(FM)/d D 

A: 0.000396489658597 
1 C: 0.000396489658593 


Thickness 


A: 0.002169495035056 
C: 0.002169495035076 


Camber 


A: 0.004203140874745 
C: 0.004203140874793 


Table 2. Figure of merit before and after single point 
designs, no thickness constraints. 


0 

Initial 

FM 

Final 

FM 

A FM 

Percent 

Change 

10° 

0.693 

0.748 

0.055 

7.9% 

12° 

0.718 

0.758 

0.040 

5.6% 

14° 

0.730 

0.761 

0.031 

4.3% 
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Table 3. Figure of merit before and after single point 
designs, thickness constraints included. 


0 

Initial 

FM 

Final 

FM 

A FM 

Percent 

Change 

10° 

0.693 

0.732 

0.039 

5.6% 

12° 

0.718 

0.747 

0.029 

4.0% 

14° 

0.730 

0.751 

0.021 

2.9% 


Table 4. Figure of merit before and after multi-point 
optimization based on the linear combination of objective 
functions given by Eq. (17). Values in parenthesis 
represent results on the refined grid. 

0 

Initial 

FM 

Final 

FM 

A FM 

Percent 

Change 

10° 

0.693 

(0.734) 

0.737 

(0.776) 

0.044 

(0.042) 

6.3% 

(5.7%) 

12° 

0.718 

(0.758) 

0.748 

(0.785) 

0.030 

(0.027) 

4.2% 

(3.6%) 

14° 

0.730 

(0.768) 

0.752 

(0.787) 

0.022 

(0.019) 

3.0% 

(2.5%) 


Table 5. Figure of merit before and after multi-point 
optimization based on the KS function given by Eq. (18). 
Values in parenthesis represent results on the refined 
grid. 

0 

Initial 

FM 

Final 

FM 

A FM 

Percent 

Change 

10° 

0.693 

(0.734) 

0.732 

(0.772) 

0.039 

(0.038) 

5.6% 

(5.2%) 

12° 

0.718 

(0.758) 

0.744 

(0.783) 

0.026 

(0.025) 

3.6% 

(3.3%) 

14° 

0.730 

(0.768) 

0.748 

(0.785) 

0.018 

(0.017) 

2.5% 

(2.2%) 


Table 6. Figure of merit before and after multi-point 
optimization based on the explicitly constrained 
approach. Values in parenthesis represent results on the 
refined grid. 


0 

Initial 

FM 

Final 

FM 

A FM 

Percent 

Change 

10° 

0.693 

(0.734) 

0.735 

(0.773) 

0.042 

(0.039) 

6.1% 

(5.3%) 

12° 

0.718 

(0.758) 

0.748 

(0.784) 

0.030 

(0.026) 

4.2% 

(3.4%) 

14° 

0.730 

(0.768) 

0.752 

(0.788) 

0.022 

(0.020) 

3.0% 

(2.6%) 
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